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Instability of the excitation spectrum of an ordered noncollinear Heisenberg antiferromagnet (AF) 
with respect to spontaneous two-magnon decays is investigated. We use a spin-1/2 AF on a trian- 
gular lattice as an example and examine the characteristic long- and short wave-length features of 
its zero-temperature spectrum within the l/S'-approximation. The kinematic conditions are shown 
to be crucial for the existence of decays and for overall properties of the spectrum. The XX Z and 
the J-J' generalizations of the model, as well as the role of higher-order corrections are discussed. 
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A quantum many-body system with nonconserved 
number of particles may have cubic vertices, which de- 
scribe interaction between one- and two-particle states. 
In crystals such anharmonicities lead to finite thermal 
conductivity by phonons In the superfluid 4 He, the 
cubic interaction between quasiparticles result in a com- 
plete wipeout of the single-particle branch at energies 
larger than twice the roton energy Q • 

In quantum magnets with collinear spin configuration, 
e.g., in AF on a bipartite lattice, cubic terms are absent 
and anharmonicities are of higher order 0]. The cubic 
terms can appear due to dipolar interactions but in 
magnetic insulators those are weak and usually can be 
neglected. It has been gradually realized that substantial 
cubic interactions should exist in the noncollinear AFs [5, 
@, 0, @| ■ Qualitatively, such cubic anharmonic terms arise 
due to coupling of the transverse (one-magnon) and the 
longitudinal (two-magnon) fluctuations in these systems. 

The noncollinearity of an antiferromagnetic spin con- 
figuration can be either induced by external magnetic 
field [|[ or by frustrating effect of the lattice {e.g., in 
the triangular lattice (TL) spins form the so-called 120° 
structure [E, 0|)- In the former case, spontaneous de- 
cays are allowed above a threshold field H* , such that 
magnons become strongly damped throughout the Bril- 
louin zone (BZ) On the other hand, the role of 

magnon interactions in the spectra of frustrated AFs is 
not well understood. The earlier work on TLAF [7( has 
discussed only renormalization of the spin-wave veloci- 
ties. The recent series expansion study [l(| has found 
a substantial deviation of the spectrum from the linear 
spin-wave theory (LSWT) and interpreted it as a sign 
of spinons. The latest work [ll| has questioned this hy- 
pothesis by showing that 1 / S expansion strongly modifies 
the LSWT spectrum leading to an overall agreement with 
the numerical data. However, the subject of spontaneous 
decays has been hardly touched upon. Astonishingly, in- 
stability of the single-particle spectrum in the presence 
of a well-defined, magnetically ordered ground state is, 



perhaps, the single most striking qualitative difference of 
the non-collinear AFs from the collinear ones. 

In this Letter, we shall study magnon decay in non- 
collinear quantum AFs at T = using an example of the 
Heisenberg AF on a TL. The noncollinearity is necessary 
but not sufficient for decays. In addition, the energy and 
momentum must be conserved within a decay process 
(kinematic conditions). Thus, the decays are determined, 
in part, by the shape of the single-particle dispersion that 
may, or may not, allow spontaneous decays. We analyze 
singularities in the two-magnon continuum that outline 
instability regions or lead to discontinuities in the single- 
particle spectrum. Our long wave-length analysis yields 
definite asymptotic statements regarding the life-time of 
magnetic excitations. We discuss briefly the XXZ and 
J-J' models where noncollinearity is also present. 

We begin by rewriting the spin-S 1 , nearest-neighbor 
Heisenberg Hamiltonian for the TLAF into the local ro- 
tating frame associated with the classical 120° structure 
of the spins, and proceed with the standard Holstein- 
Primakoff transformation of the spin operators into 
bosons followed by the Bogolyubov transformation diag- 
onalizing the harmonic part of the bosonic Hamiltonian. 
This procedure leads to the following Hamiltonian: 
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where the ellipses stand for the classical energy, other 3- 
boson terms that do not lead to decays, 4-boson, and the 
higher-order terms. Although we will need some of these 
other terms for the 1 /^-expansion below, Eq. |TJ) will 
suffice for the purpose of this paper. All of the necessary 
terms can be found in 0, 0] • The LSWT magnon energy 
and the 3-boson vertex in Eq. (jTJ) are given by: 



£k = 3JS\/(l- 7 k)(l + 2 7k ) , 

Vk,q = 3l J \/'SS/2 [/q,q',k + /q',q,k " 

where J is the exchange constant, q' = k 
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FIG. 1: (color online). Brillouin zone of the TL. The 
shaded area corresponds to the region where spontaneous 
two-magnon decays are allowed. T = (0,0), K=(47r/3, 0), 
K'=(27r/3,27r/ v / 3), Y=(0,ir/y/S), and M=(tt, tt/V3) points 
are highlighted. The lines correspond to the extrema in the 
two-magnon continuum described in the text. 

7k = l(cosk x + 2cos^cos^§^), % = l(sink x - 

2 sin 4f cos ) , and u, v are the Bogolyubov coeffi- 
cients: 2uf - 1 = 3JS(l + 57f)/ei, uf - t>? = 1. 

Kinematics: long wave-length limit. — The magnon 
branch in Eq. |2|) has three zero-energy modes: k = 
and k = ±Q= (±§tt,0), points T, and K (K') in Fig. [TJ 
respectively. In contrast with a square-lattice AF, veloc- 
ities of these Goldstone modes are different: cq — \/3/4 
and c±q = \/3/8 (in units of 3JSa). This immediately 
implies that excitations with k — > are kinematically 
unstable towards decays into (q, q') — > (Q, — Q) ones, in 
close analogy with a decay of a longitudinal phonon into 
two transverse ones [l[ . Clearly, such decays are immune 
to 1/5 corrections as long as velocities remain different. 
This picture is pertinent to all other non-collinear AF 
with more than one Goldstone mode. 

For magnons at small k = k — Q there exists a more 
subtle reason for decays. Instead of the usual convex 
and isotropic form, the magnon energy is nonanalytic 
with varying convexity: ~ cgk(l — a v k), where 
a v ~ cos 3tp. This form together with the commensura- 
bility of the ordering vector create kinematic conditions 
for the decays from the steeper side of the energy cone at 
k — * Q into the less steeper sides at q, q' — * — Q. Thus, 
magnons near the Q-point are unstable only in a range 
of angles. Although such conditions are more delicate, it 
is very unlikely that the higher-order terms would selec- 
tively cancel the nonanalyticity. Therefore, the k — > Q 
decays should be prominent in the TLAF. 

Kinematics: full BZ. — In the model (JTJ an excitation 
with the momentum k is unstable if the minimum energy 
of the two-particle continuum, -Ek.q = £q + £k-q, is lower 
than £k- Then the boundary between stable and unstable 
excitations is where such a minimum crosses the single- 
particle branch and the decay condition -Ek,q = £k is first 
met. Thus, to find these boundaries one should analyze 



the extrema of the continuum. For a gapless spectrum 
there can be several solutions as we show for the TLAF: 

(a) Decay with emission of a q = magnon. 2?k,o = £k 
for any k but it never crosses the magnon branch. 

(b) Decay with emission of q = ±Q magnons. Equa- 
tion that defines the boundary is E^q = £k±Q = £k and 
its solution is shown by the dotted line in Fig. [TJ The 
shaded area is where magnon decays are allowed. It can 
be shown that i?k,Q corresponds to an absolute minimum 
of the continuum within the shaded area. In accord with 
our long wave-length discussion, the area around k = 
is enclosed and it is a finite segment in the vicinity of the 
Q-point where decays are allowed. 

(c) Decay into two identical magnons. The two-magnon 
continuum has extrema that are found from dE^^/dq = 
@, @| • This condition means that the products of de- 
cay have equal velocities. The simplest way to satisfy 
that is to assume that their momenta are also equal. 
This is fulfilled automatically if q = (k + Gi)/2, where 
Gj is one of the two reciprocal lattice vectors of a TL, 
Gi = (2 7 r,27r/ v / 3) and G 2 = (0,4tt/V3). The curves 
for the solution of ek = 2e( k +Gi)/2 are shown in Fig. [JJ 
by the solid lines. Note, that in the case of the TLAF 
these extrema are not the minima, but the saddle-points. 
Nevertheless, they lead to essential singularities in the 
spectrum as will be discussed below. 

(d) Decay into non-identical magnons. In a more general 
situation, the decay products with the same velocities 
may have different momenta and energies £k-q =/= £q- 
Then, one has to solve the decay condition i?k,q = ^k to- 
gether with the extremum condition dE-^^/dq = 0. The 
solution is shown in Fig. [JJby the dashed line. As in (c), 
corresponding extrema are the saddle-points. 

Thus, the area of two-magnon decays in Fig.[T]is deter- 
mined by the solution (b) as it encloses regions (c) and 
(d). This may not be the case for other systems (see be- 
low the XX Z model). Generally, the area of the decays 
is a union of the regions given by (b), (c), and (d). 

Spectrum: 1/5 corrections. — The 1/5-correction to 
the TLAF magnon spectrum is given by the one-loop 
self-energy diagrams from the 3-boson terms, and the u>- 
independent 4-boson contribution, &k = £k(£k) + Se^: 

V f ^ i^/ ^ql 2 l^ql 2 ^ 

2 q \^k( w ) ^kM / 
fcf = 9 J 2 S [Anl + A 27k - A 3 ] /4e k (5) 

where D^(uj) = ui =p e q =p £k-q ± £0, the "source" 3-boson 
vertex is V" k ,q = 3iJy/3S/2 [g q ,q',k + ffq'^k - 3k,q,q'], 
and Ai = (4c + c x - 5c 2 - 4), A 2 = (~2c a + ci+ c 2 + 2), 
A3 = (2co+ 2c\ — 4c 2 — 2). Here, wc have defined the 
following constants: c„ = 3 J5 X]k 7k/ £ k, Co = 1. 574733, 
ci = -O.f 04254, c 2 = 0.344446. 

We calculate the spectrum £k = £k + fek for 5 = 1/2 
using numerical integration in Eq. (J^EJ) . Fig.[5]shows £k for 
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FIG. 2: (color online). Upper row: magnon dispersion along 
the lines TK and YM in the BZ, see Fig. [Tj The dashed lines 
are the LSWT spectrum, the dotted lines is the bottom of 
the LSWT two-magnon continuum, and the solid lines are 
the spectrum with the 1/S correction. Lower row: imagi- 
nary part of the 1/S magnon energy along the same lines, kb 
is the intersection point of the one-magnon branch with the 
two-magnon continuum along the YM line, k* and k! points 
correspond to the singularities discussed in text. 



two representative directions in the BZ together with the 
LSWT spectrum and the bottom of the two-magnon con- 
tinuum. The l/jS-renormalization of the magnon spec- 



trum is quite substantial, see also [11] . This is because 
the continuum strongly overlaps with, or, for the k-areas 
outside the decay region, has significant weight in a close 
vicinity of, the magnon branch. Such a purely kinematic 
effect explains a mystifying dichotomy: quantum correc- 
tions to the spectrum in the TLAF are large compared 
to the square lattice AF, while the ordered moments are 
about the same [H, 0, E3] • 

Decays: long wave-length limit. — The decay vertex 
© for magnons near the T point scales as Vk,Q+q oc 
(<?' — q)y/k/qq', for small q and q' = |k — q|. A simple 
power counting yields the leading term in the imaginary 
part of the self-energy. In a typical decay q,q' ^k giving 
for the decay probability: | Vk,Q+q| 2 cx k. Since there is no 
constraint on the angle between k and q, the 2D phase 
volume restricted by the energy conservation contributes 
another factor of k, such that 9m{Ek(£k)} oc k 2 . A more 
detailed analytical calculation yields Sm{Ek(ek)} ~ 
— 0.789 Jk 2 , in agreement with the data in Fig. [2] 

The decay vertex for k — * Q magnon has a more con- 
ventional scaling: Vcj+k.-Q+q oc ^Jkqq' , so the decay 
probability is lypocfc 3 . Due to a constraint on the angle 
between k and q, the decay surface in q-space is a cigar- 
shaped ellipse with length ~ k and width ~ fc 3 / 2 that 
makes the restricted phase volume of decays to scale as 
fc 1 / 2 . This results in a nontrivial fc 7 / 2 scaling of the de- 



cay rate. Numerically, along the TK line 3m{Ek(£k)} ~ 
— 1.2J&; 7 / 2 . In a similar manner, one can show that at 
the boundary of decay region (e.g., point k b in Fig. [5]) 
the decay rate grows as 9m{E k (£k)} oc (fc — k b ) 2 . 

Spectrum: singularities due to decays. — A remarkable 
feature of the spectrum in Fig. [2] is the singularities in 
the real and the imaginary parts of £k- Clearly, they are 
due to spontaneous decays as it is only the decay term 
in ([4]) that contributes to the imaginary part. 

A close inspection shows that k' and k* singularity 
points in Fig. [2] correspond exactly to the intersection of 
TK and YM lines with the saddle-points in the contin- 
uum, solid and dashed lines in Fig. [TJ The "strong" (k*) 
and "weak" (fc') singularities correspond to decays into 
identical and non-identical magnons, solutions (c) and 
(d) above, respectively. Fig. [3] shows the decay contours, 
i.e., ID surfaces in a 2D q-space into which a magnon 
with the momentum k can decay, for kl and k* saddle 
points along the TK line. In both cases, these contours 
undergo topological transition. 

Close to such a transition the denominator in Eq. (j4|) 
is expanded as £k— -Ek, q ~(th— v 2 )^k— (3 x q 2 +fl v q 2 , v\ and 
V2 are velocities of the initial and final magnons, respec- 
tively, j3i are constants, Ak = k—k*, k* is the saddle point. 
Integration in Eq. ((4|) yields a logarithmic singularity in 
the imaginary part SmjEi} oc — ln|A/Afe| and a con- 
comitant finite jump in the real part of the self-energy 
5Re{Efe} oc sign(Afc), A is a cutoff. The cutoff (size of 
the "bubble" in Fig. [3]) is small in the case of "weak" 
singularities, which explains the weakness of them. The 
logarithmic form agrees perfectly with the data in Fig. [2] 
Overall, our analysis of the long- and short wave-length 
behavior gives complete understanding of the 1/5 results 
for the magnon spectrum in TLAF. 

An important question is whether the singularities in 
the spectrum will withstand the higher-order treatment. 
The first possibility is when at least one of the final 
magnons is itself unstable. Then the log-singularity will 
be cut off by the decay rate of the final magnon and 
QmjEfc} will have, at most, a weak maximum near the 




FIG. 3: (color online). The upper right parts of the BZ in 
the q-space. The decay contours for k = k' (left) and k — k* 
(right) along the TK line are shown. Both k-points correspond 
to the saddle points in the two-magnon continuum. The cor- 
responding q-contours undergo topological transition. 
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FIG. 4: (color online). Decay regions and singularity lines 
for the XXZ model, a — 0.96 (left) and the J — J' model, 
J'/ J — 0.9 (right). Definition of lines is the same as in Fig.[T] 



topological transition. For the TLAF such a scenario 
is realized for a large fraction of the "weak" singular- 
ities (k' in Fig. [2]). However, all of the "strong" sin- 
gularities (k* in Fig. [2]) and some of the "weak" ones 
belong to another class, in which both magnons created 
in the decay are stable. We have checked that the first- 
order l/S corrections do not shift appreciably the insta- 
bility boundaries in Fig. [1] Hence, at the saddle points, 
the logarithmic divergence of the one-loop diagrams will 
persist even for the renormalized spectrum. In such a 
case, vertex corrections become important 0. Summa- 
tion of an infinite series of loop diagrams yields the self- 
energy from the decay processes near the singular point: 
£k(w) — ia/\n\Auj — i^Afe], with a > 0. One may con- 
clude that the decay rate becomes vanishingly small as 
Aw, Ak — > 0. This is, however, not true. An attempt 
to solve the Dyson's equation Gj^ 1 (o-') = with this self- 
energy yields no solution for uj near the real axis. There- 
fore, the decay rate of quasiparticles around solid lines in 
Fig. [T] will remain large and quasiparticle peaks will be 
strongly suppressed even for large S. 

Other models on a TL. — Two straightforward gener- 
alizations of the Heisenberg model on a TL are (i) the 
anisotropic XXZ model with a = J z /J xy < 1 and (ii) 
the J-J' model for an orthorhombically distorted tri- 
angular lattice with one of the interactions within the 
triangle (J) stronger than the other two (J'). 

In the spectrum of the XXZ model magnons at Q are 
gapped with the gap £q cx a/1 — ol at a < 1. This has 
two immediate consequences: (a) magnons at k — > are 
stable up to £k = 2eq and (b) Q-magnons become stable 
themselves. Thus, the star-shaped decay region in Fig. Q] 
develops a hole in the middle and have vertices shrunk 
and rounded. The evolution of the decay boundary with 
a < 1 is non-trivial. Initially, the emission of a Q magnon 
remains an absolute minimum of the two-magnon contin- 
uum for most of the decay region. At a± « 0.993 the decay 
into non-equivalent magnons switches from being a line 
of saddle points into the absolute minima of the contin- 
uum and takes over the decay boundary. Fig. [¥] shows 
the instability and the singularity lines for a representa- 
tive value a — 0.96. Further decrease of a completely 



eliminates the decay region at around a w 0.92. Thus, 
magnon decays are present in an anisotropic TLAF, but 
only at not too large anisotropies. 

For the J-J' model the Goldstone modes at k = ±Q 
are preserved but the ordering wave-vector Q becomes 
incommensurate. This does not change the kinematics 
of the decays for the k — > magnons, but forbids the 
decays from the vicinity of the Q point into the vicinity of 
— Q point as the quasi- momentum cannot be conserved. 
However, the decays in the vicinity of (inequivalent now) 
K' points are still allowed. Overall, the decay region 
grows with the decrease of J'. At J' » 0.34J, relevant 
to CS2CUCI4 [12J, the decay region covers most of the 
BZ. With the decrease of J' the LSWT single-magnon 
dispersion develops a low-energy branch in the direction 
perpendicular to the "strong" J. That makes the rest of 
the spectrum prone to decays into it. 

Conclusions. - We have shown that magnon de- 
cays must be prominent in a wide class of noncollinear 
AFs. We calculated the decay rate in the spin-1/2 TLAF 
within the spin-wave theory. In the long-wavelength 
limit, the life-time of low-energy excitations is predicted 
to exhibit a non-trivial scaling. For the short-wavelength 
magnons, the decay rate is large, 29m{ek} ~ 0.45fte{ek}, 
in a substantial part of the BZ. Topological transitions 
of the decay surface also lead to strong singularities in 
the spectrum that remain essential even for large values 
of spin. Therefore, excitations in ordered, spin-5 1 , AFs 
may not necessarily be well-defined for all wave-vectors. 
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